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Abstract 

We have implemented the so called GW approximation (GWA) based on an 
all-electron full-potential Projector Augmented Wave (PAW) method. For the 
screening of the Coulomb interaction W we tested three different plasmon-pole 
dielectric function models, and showed that the accuracy of the quasiparticle 
energies is not sensitive to the the details of these models. We have then 
applied this new method to compute the quasiparticle band structure of some 
small, medium and large-band-gap semiconductors: Si, GaAs, AlAs, InP, 
SiMg2, C and (insulator) LiCl. A special attention was devoted to the con- 
vergence of the self-energy with respect to both the k-points in the Brillouin 
zone and to the number of reciprocal space G-vectors. The most important 
result is that although the all-electron GWA improves considerably the quasi- 
particle band structure of semiconductors, it does not always provide the 
correct energy band gaps as originally claimed by GWA pseudopotential type 
of calculations. We argue that the decoupling between the valence and core 
electrons is a problem, and is some what hidden in a pseudopotential type of 
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I. INTRODUCTION 



In the last few years, with the increase of computer power, many researchers have tried 
to unify two major but different approaches to the computation of electronic properties of 
materials: (1) The local density approximation (LDA) to the Kohn-Sham equationsS where 
the correlation effects are included in an average way through the use of the parameterized 
local exchange-correlation potential.! (2) Many-body models based on simplified and param- 
eterized Hamiltonians omitting much of the subtlety of the chemical bonding but handling 
the correlation effects effectively. The aim of the unification being to provide a good de- 
scription of the quasi-particles properties in real solids without omitting either the subtlety 
of chemical bonding or the correlated nature of the electrons.iHHil 

That LDA is inadequate has been known for example in the case of semiconductors: 
the LDA band gaps are at least 50% smaller than experimental values. It was shown by 
Perdew and Levyl! and Sham and SchliiteiH that the difference between the highest occupied 
and lowest unoccupied state of the density functional theory (DFT) eigenvalues is not the 
true quasiparticle band gap, but differs from it by A, the discontinuity in the exchange- 
correlation potential when an electron is added to the system. Later, it was shown that 
quasiparticle calculation based on the the GW approximation (GWA) of Hedin&i produced 
band gaps which were in good agreement with experiment .EHll This development in the 
theoretical study of the electronic structure of solids is the inclusion of many-body effects 
in the calculation, principally, through the computation of the quasiparticle self-energy and 
the results are an adjustment of the energy splittings obtained within the local density 
approximat ion .0 

After this initial success the GWA has been extensively used to study other properties 
of different type of materials: (1) The band-width narrowing in alkali-metals, and their 
clusters. (2) The surfaces states of semiconductors, i.e., improving band gaps of surface 
states,!! looking for dimer buckling in Si surfaces.!! (3) The effects of correlation on the 
valence off-set between different bulk semiconductors.il (4) The character of the band gaps 
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of superlattices and anisotropy of optical matrix elements.El (5) The orientational disorder 
and photoemission spectra of solid (6) The electronic properties of elemental Ni and 

its energy loss spectra.0 (7) The quasiparticle properties of atoms using various GWA.i (8) 
The Schottky barrier between a metal and a semiconductor.0 (9) The inclusion of excitonic 
effects in the calculation of the dielectric function of semiconductorsil~0. All these different 
studies in a very short time established the GWA as a good 'first-principles' method for 
computations of quasiparticle properties of real materials. 

In most of these studies the plane wave expansion of the LDA pseudopotential basis- 
set is used, which makes an analysis in terms of chemically relevant orbitals difficult. An 
additional difficulty with doing the GWA to the self-energy in a plane wave basis set is 
that the computational effort for systems with localized electrons is enormous since the 
plane wave expansion of the wave function becomes a real problem especially for systems of 
localized electrons or low dimensional materials. 

In this paper, we propose a new implementation of the GWA based on an all- 
electron method using the recently developed all-electron projector augmented wave (PAW) 
method.il Our constructed ab-initio method of quasiparticles is then based on the framework 
of the LDA theory in conjunction with the PAW method. The knowledge of the one-electron 
Green function provided by the PAW method allows us to construct the quasiparticle self- 
energy within the GWA, in which the dynamical screening of the electron-electron interaction 
arises from a plasmon model dielectric functionS'ElS for which the parameters are adjusted 
to the dielectric function calculated in the random- phase-approximation0 (RPA). The full 
calculation of systems of many atoms in the unit cell is prohibitive on sequential comput- 
ers, so we had to parallelize our numerical code so that the calculations become feasible on 
scalable parallel architectures. 

Our paper is organized as follows: In the fist part we introduce the projector augmented 
wave method which is used to solve the Kohn-Sham equations and provides the one-electron 
Green's function and the RPA dielectric function which are then used to compute the self- 
energy. We then describe our GW implementation based on the PAW method and show 



the difficulties, traps, encountered when implementing an all-electron GWA. In particular, 
we mention the difficulty related to the decoupling of the core and valence electrons. We 
discuss also in details the use of symmetry to reduce the computational time of the self- 
energy and the dielectric function, as well as the parallelization of our numerical code. In 
the third section we apply our method to determine the electronic structure of two distinct 
semiconductor groups: some small and medium-band-gap semiconductors: Si, GaAs, AlAs, 
InP, Mg2Si and some large-band-gap semiconductors (insulator): C, and (LiCl). We then 
compare our results with available GWA calculations and experiments. 



II. METHOD OF CALCULATION 

A. A brief introduction of the PAW method 

In the density functional theory implemented in the framework of the LDA0, an electronic 
structure calculation requires the solution of Kohn-Sham type of equations in a self consistent 
way. The computation of material's band structure consists in finding the Bloch wave 
functions ty n k(r), where n and k denote a band index and a wave vector in the Brillouin zone, 
respectively. In the projector augmented plane wave (PAW) formalism,il all calculations 
are performed on a smooth "pseudo" wave function \l/„k(r) which is expressed as a linear 
combination of plane waves. The passage from the smooth pseudo-wave function to the 
all-electron wave function exhibiting the correct nodal behavior in the augmentation regions 
(spheres centered on each atom) is achieved by defining three atomic type of functions in 
each augmentation region: (1) The all-electron basis functions $"(r), (2) The 'pseudo' basis 
functions $"(r), (3) The projector functions p"(r). Here i = k, rrii, rii where U and rrii denote 
the orbital and magnetic quantum numbers, respectively. The index rii is introduced to leave 
the opportunity of choosing more than one function for each angular momentum channel 
(lj,mj). These functions are defined so that: 

$°( r ) = $«( r ) for r > r a c (1) 
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where r a c are the radii of non-overlapping spheres centered at each atomic site a. The 
projector functions vanish for r > r a c and satisfy the orthogonality property: 

m*°) = hi (2) 

Using these functions, the all-electron wave function ^/ n k(r) can be deduced from the 'pseudo' 
wave function ^ nk (r) according to the relation \l/ n k( r ) = ^rik( r ) + ^nk( r ) ~ ^n.k( r ) with 

*nk(r) - *ik(r) = £[ - * a ) - *?(r - R a ) ]<p?|*„k> (3) 

a,i 

where R a denotes the atomic position of the atom a in the unit cell. It is useful to 
point out the fact that \l/^ k — ^ k vanishes in the interstitial region and defines the quantity 
necessary to describe the true wave function in the augmentation regions while \I/ n k describes 
the true wave function in the interstitial region. 

The PAW formalism is designed to easily calculate the expectation value of local or semi- 
local observables. For example the expectation value of an operator A(r) between two Bloch 
wave functions \I/ n k and \l/ m k can be calculated as a sum of three contributions: 

Aik.mk = (^nkl^l^mk) = ^nk,mk + A^ k ,mk — ^nk,mk (4) 

where 

ink,mk = (^nkl^l^mk) (5) 

This contribution is evaluated in the plane wave basis set. The last two contributions 

= E<*»*I^>< $ ?I^I*?><P?I*»*> ( 6 ) 

i,j,a 

and 

= Y,^n k mmA\^)(p a ^ mk ) (7) 

are evaluated in the augmentation regions. 

Blochli managed to decompose the total valence energy into three contributions in complete 
analogy with the decomposition of the expectation value of an operator 
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E = E + E 1 -E 1 (8) 

where E includes the kinetic energy of the valence electrons, the interaction energy with 
the nuclei having atomic number Z a and with the core electrons, which are described in the 
frozen core approximation, the Hartree energy and the exchange correlation energy of the 
valence electrons. 

Minimizing the total energy with respect to the (variationnal principle) leads to a 
generalized eigenvalue problem which is solved in a self-consistent way, giving the pseudo 
wave functions from which the all-electron wave functions \l/„k are easily deduced by 
means of equation |3|. 



B. The GW approximation 

1. Quasiparticle calculation in the GW approximation 

Solving the Kohn-Sham equations, where the exchange-correlation effects are included 
in a mean way in an exchange-correlation potential V^P (r) yields to eigenvalues which 
can not be assimilated to the excitation energies of the solids. Indeed, there isn't any 
Koopman's theorem for the Kohn-Sham theory and these eigenvalues only have meaning as 
Lagrange parameters. The Green's function theory is a good tool for properly describing the 
excitation energies. In the quasiparticle approximation, we can find the excitation energies 
of the system by solving a quasiparticle equation 

(T + + Wkn(r) + J d 3 r'£(r, r', E n (k))fa n (r') = E n (k)^ n (r) (9) 

instead of locating the poles of the Green's function. Here, T is the kinetic energy operator 
(— |V 2 in atomic units), V ext is the external (ionic) potential, Vh is the Hartree potential due 
to the average Coulomb repulsion of the electrons and £ is the self-energy operator which 
summarizes the many-body effects. £ is in general a non local, energy- dependent, non- 
Hermitian operator. The non-Hermitian part of E gives rise to complex eigenvalues E n (k). 
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The real part of the eigenvalue is associated with the energy of the quasiparticle while the 
imaginary part is related to the inverse of the lifetime of the quasiparticle. Fortunately, 
the imaginary part of the eigenvalue, which can be attributed to the possibility for the 
quasiparticle to lower its energy via particle-hole formations (the decay of the quasiparticle 
via emission or absorption of plasmons is not effective near the Fermi level), becomes small 
near the Fermi level because of the restriction of the phase space for these processes. This 
is known as Pauli blocking. So, the quasiparticle concept is effective near the Fermi Level. 

The similitude between the quasiparticle equation and the familiar self-consistent field 
equation in the Kohn-Sham formulation of the density-functional theory 



is striking if we set J] LDA = 5(E)5(t,t')V xc (t) in the quasiparticle equation. So the Kohn- 
Sham equation can be regarded as an approximation to the quasiparticle equation where 
the self-energy is approximated by a local and energy-independent potential. 

A more realistic, but relatively simple approximation to the self-energy, which takes 
into account both non-locality and dynamic correlations, known as the GW approximation 
(GWA), was initiated by HedirJU. This approximation was originally derived from a many- 
body perturbation theory as the first term in the expansion of the self-energy in the screened 
interaction W. Within this scheme, the self-energy 



is approximated by a convolution with respect to the frequency variable of the Green's 
function with the screened interaction along the real axis. Here 5 is a positive infinitesimal 
which ensures that only the poles of the Green's function associated to occupied states 
contribute to the integral, G is the best available Green's function and W is the best available 
screened Coulomb interaction. Accordingly, G is taken to be the Green's function built from 
LDA orbitals 




(10) 




(11) 



= lim 

nk 




(12) 



e n (k) + i5sgn(e n (k) - ji) 



8 



where \x denotes the chemical potential, and where 5 — > + means 5 being a positive in- 
finitesimal number. 

The dynamically screened interaction is defined by 

W{r,r',u) = /"dr'V^forVMrVO (13) 

where v denotes the bare Coulomb interaction and e the dielectric function defined by 

e(r, r', oj) = S(r - r') - / v(r, r")P°(r", r', u)dr" (14) 

where the irreducible polarisability P° is calculated within the RPA approximation (Bubble 
approximation) 

P°(r, r', u) = ~ f du'Gir, r', co)G(r', r, J - u)e i5 "' (15) 

Z7T J 

Using the LDA expression of the Green's function and integrating in the complex plane leads 
to the familiar expression of the irreducible polarisability 

p o (r _, x = V MO - n m (kO]vl> k; (r)^(rOvl> k , m (rOvl>^ m (rO 

1 ' ' ] e,(k) - e m (k') -u + i6[m(k) - n m (k')} { ] 

After having defined the different quantities, which enter the self-energy in the GW approxi- 
mation, it is natural to look for the solution of the quasiparticle equation along the following 
line 

^km( r ) = J2 a ™n(k)^ kn (r) (17) 

n 

It has been shown by inspection that, a mn ~ 8 mn , which means that the quasiparticle 
Hamiltonian is virtually diagonal in the \l/kn basis for semiconductors such as Si!. A plausible 
argument for neglecting the non-diagonal part has been given recently by Hedin.S Since the 
quasiparticle wave functions ^ are similar to the LDA wave functions ^/^ n (the overlap 
between the quasiparticle wave function and the LDA wave function is about 99 %i), the 
numerical work is therefore considerably reduced. Indeed, using first order perturbation 
theory yields the following result: 



E n (k) = e„(k) + (^ kn |£(r,r',£ n (k))|^ kn > - (^ kn |V; L c DA (r)|^ kn ) (18) 

In principle, the solution to this equation should be obtained via an iterative method, but 
expanding Eq. (|18D to first order in energy around e n (k) yields accurate numerical results. 
Then the quasiparticle energy can be obtained via 

E n (k) - e n (k) = Z nk [(* kn |S(r,r',e n (k))|* kn ) - (* kn |F^(r)|* kn )] (19) 

where the quasiparticle renormalization factor Z nk is 

Z nk = [1 - (^^(r^e^k))^)]- 1 (20) 

Since £S < 0, we have < Z„ k < 1 with typical values of 0.8 for bands close to the band 
gap and for all materials considered here (see Table |). Values of Z nk about 0.8 imply that 
we still have well-defined quasiparticles in the system but that 20 percents of the spectral 
weight is now distributed over a range of frequencies. 

2. Calculating the matrix elements of £ within the PA W formalism 

As can be seen from Eq. (|19"D, the central problem of this scheme consists in evaluating 
the diagonal matrix elements of the self-energy between LDA orbitals. The quantities which 
enter the self-energy are functions / of two locations r and r'. These functions have the 
following translational symmetry property f(r + R, r' + R) = f(r, r') where R is a Bravais 
lattice vector. Then, we can fix the Fourier transform convention for such functions 

/(r,r» = l J2 e^ +G )-V GG '(q^)e- j(q+G ' ) - r ' (21) 

^ q,G,G' 

where q is a wave vector in the first Brillouin zone, G a reciprocal lattice vector and Q the 
crystal volume. 

In the plane wave basis, the dielectric function can be defined as 

ecc (q, oj) = 8 GG > - u(q + G)P GG , (q, u) (22) 
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Then according to Eqs. ([□] and [L3D, the dielectric matrices eGG'(q, have to be calculated 
and inverted for many values of u. This is computationally time consuming. Nevertheless, 
it has been carried out by some authorsi^ who choose to evaluate the frequency integral by 
using a Gaussian integration scheme along the imaginary axis to circumvent the problem of 
the pole structure of the screened interaction along the real frequency axisi. An alternative 
approach is the use of a plasmon pole modefi0'0S to mimic the frequency dependence of 
the dielectric matrix. These models give a good description of the low frequency behavior of 
the dynamically screened interaction and allow the determination of an analytic expression 
for the frequency integral appearing in Eq. fllTl) . We used three types of plasmon pole 
models to describe approximately the dependence of e _1 (u;) on the frequency uj. These 
different models were proposed by Von der Linden and Horsch0, Engel and FaricS, and 
by Hamada et alM We choose to detail the first one and emphasize that the values of the 
quasiparticle energies are not sensitive to the model type (see Table [Ij which details the 
calculated QP energies of Si with these three models). 

In such an approach, we consider the symmetrized dielectric matrix eGG'( ( l> a; ) defined 

by 

cgg' (q, w) = r^TTTi ecc (q, w) (23) 
|q + yjr | 

Introducing the following notation 

M£ m (k, q) = (v& k _ qn |e-^+ G >- r |vI> km > (24) 



and using Eqs. (|23| , [22j , |16|) leads the following result for the static symmetrized dielectric 
matrix 

~ , 167T M G -(k,q)[M G c ,(k,q)]* 

(C1 ' U ~ 0) " 6GG ' ' n|q + G||q+G'| J k e ,(k-q)- ec (k) (25) 

We then diagonalize the hermitian e(q, u = 0). Its real eigenvalues A p (q) and orthonormal 

eigenvectors |0 p (q)) can be used to perform the matrix inversion according to 

r\ q ,uj = 0) = £ |0 P (q)) Y | T (0 P (q)| (26) 

n 



Von der Linden and Horsch supposed that the frequency dependency of the dielectric matrix 
is encompassed in the eigenvalues of the static symmetrized dielectric matrix while the 
eigenvectors are frequency independent 

^G'(q^) = EC(?)T 7 L ^[^(q)r (27) 



where 



Ap(q,a;) 1+ 2 ^ p(q){ uo - (u p (q) - i5) to + (u p (q) - i5) } (28) 
It should be noted that such a parameterization of the frequency dependence of the inverse 
of the eigenvalues guarantees that e is an even function of u. Here, two parameters have to 
be determined: (1) The force of the pole z p (q) which is defined by letting u = in Eq. (|28|) 

^(q) = l-A p 1 (q,0), (29) 



and is positive since it can be shown that A„ 1 (q, 0) lies in the interval (0, 1)EH (2) The 
frequency of the pole a> p (q) which is determined from the Johnson's sum ru lei. If we 
introduce the following quantities 

e G (q) = (30) 

and 

L GG , = (q + G)(q + G>(G - G'), (31) 

where p(G) denotes a Fourier component of the charge density, we can write the frequencies 
of the poles according to 

Air 

^? = -^mci)\L(q)m q )) (32) 

Using the usual development of p and 9 along the reciprocal vectors, it can be easily shown 
that 

(e*(q)|L(q)|e*(q)> = \\ p(r)|V6(r)| 2 rf 3 r (33) 
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where V denote the volume of the unit cell. Since z p (q) and the quantity defined by Eq. 
(|33|) are both positive we deduce that the frequency defined by Eq. fl32|) is positive. Now 
all the ingredients of the plasmon pole model have been defined, and we are able to give 
an analytic expression of the matrix elements of the self-energy in the GW approximation 
between LDA Bloch orbitals. 

To test our implementation of the plasmon pole model, we have plotted in Fig. p] the 
Engel-Farid Plasmon model band structure of Si along L, T, and X high symmetry directions. 
We have found that our implementation is in excellent agreement with the results of Engel 
and Faridil, and Aulbur0. Indeed, for small k wave vectors, the lowest plasmon band shows 
a quadratic dispersion u; (k) = u (0) + «|k| 2 , with a dimensionless direction-dependent 
dispersion coefficient a (see up-triangle curve in Fig. |l|). We find u (0) =15.7 eV and 
«a =0.33 in good agreement with the values of 15.91 eV and 0.34 of Engel and Farid as well 
as the experimental values of 16.7 eV and 0.410. To check also the validity of the plasmon 
pole model to be used as a substitute for the dielectric function of real materials we have 
compared the plasmon pole model of Hamada et al.0 with our direct ab-initio computation 
of the dielectric function within the RPA including the so called local-field effects (see Fig. 
0) and with available experimental results^ for the energy loss spectrum. We notice that 
the model mimics nicely our ab initio calculated dielectric function. 

The matrix elements of the self-energy could be divided into an energy independent 
contribution Yl 1 ? and an energy dependent contribution £(c<j). The first contribution corre- 
sponds to the Hartree-Fock contribution and is given by 

( , u ^ = _£ EEE !*jM (34 ) 

where the summation runs only over the occupied states. The energy dependent contribution 
can be expressed as 

Wta M w WW-7T 71 7T1 — ~T~\ ^ 7 71 l k ' 

n w - e ™( k - q) + K(q) - ™\sgn{ii - e m (k - q)) p 

(35) 
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where 



/?™(k, q) = Et M G n (k, q)]*e pq (G) (36) 

G 

and 



sgn(fj, - e m (k - q)) 



1 for e m (k - q) < /i 
-1 for e m (k - q) > \i 

Here m denotes an electronic band index, p a plasmon band index and q a vector in the 
Brillouin zone. It should be emphasized that the summation over m is not restricted to 
occupied states as in the expression of the Hartree-Fock contribution. Both the poles of the 
Green's function and of the screened interaction contribute to this expression. 

3. Numerical details 

One of the central problem within the realization of the GW approximation is the cal- 



culation of matrix elements whose type is defined by Eq. (24]). Using the PAW formalism, 
the smooth 'pseudo wave function' \I/ kri associated to an 'all-electron' LDA wave function 
^kn can be written 

^n = ^=Y. A UG)^ +G)r (37) 
v V 



G 

where the sum runs over reciprocal lattice vectors. As illustrated in the part dedicated to the 
PAW formalism, the expectation value of e~ l ( q+G )" r can be divided into three parts according 
to M£™(k, q) = (^ k - qn |e-^+ G )- r |^ km > + (* k _ qn | e - i(q+G) - r |* km > - ^i-^^^iJ ■ 
The first term, which involves plane waves is defined as 

£A*_ qn (G')A km (G + G') (38) 

G' 

In our scheme, the summation over reciprocal lattice vectors in this expression takes account 
of a vector G' if both G' and G + G' are smaller than a cutoff parameter. In general 300 G' 
vectors were included in the summation for the systems studied here, except for Diamond 
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where the convergence is achieved only when 400 G' were used. The two remaining terms 
which involve localized contributions can be expressed as 

£<*k-*»|P?>[ (<^|e-^+ G )- r |$f> - (<£»|e-^+ G )- r |<£»> M|* mk > (39) 

Since the overlap between the 'pseudo wave functions' and the projectors is known, we have 
to calculate quantities like ($ a |e~^ q+G )- r |$f ). Using the development of plane waves along 
spherical Bessel functions ji, we get 

•/ n, I 47r e^( q+G )- Ra Elm H)' Y lm (q + ~G) G\ m m ,._ 

($ a |e ( q+ )' r |$ a ) = < 1 (40) 

[ x f^dr r 2 j,(|q+G||r|) $/, n ,(r) $^(r) 
where the Gaunt coefficients are defined to be 



G\T mi hm] = V^JcM Y* mi (?) Y lm (r) Y ljmj (?) (41) 



The symmetrized dielectric matrix defined by Eq. fl25D as well as the matrix elements 
of the self-energy defined by Eqs. (|34|) and ([35]) are obtained by an integration over the 
Brillouin zone. We used the special points' technique^! in which the summation over a 
uniform mesh of k-points is reduced by symmetry to a summation over fewer special It- 
points if the integrand possesses the full symmetry of the point group of the lattice. The 
summation over q-points in the expectation value of the self-energy has to be carried out 
carefully, since the integrands have an integrable singularity in for q — > 0. This can be 
readily seen in the expectation value of defined by Eq. ( p4|) where the divergence occurs 
when G = 0. We follow the idea of Gygi and Baldereschi@ introduced for fee lattice by 
considering a smooth function F(q) which reflects the translational symmetry of the Bravais 
lattice and which diverges as ^ as g vanishes. If we have to integrate a function g(q) which 
behaves as 4r for small q, we can write 

E *(q) = E b(q) - A H<d) + a E ^(q) (42) 

q q q 

Such a decomposition is illustrated for the Hartree-Fock contribution 
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q m occ 



!Mr(k,g)l 2 

^ 'q+G|* " F(q)5 " 



47T 

tt E tnEf(q) (43) 

" m occ g 



G 

The function in square brackets does not contain any divergence owing to the following 
property 

limM m "(k,q) = 5 mn (44) 

q^O 

and is easily integrated by special points method while the integral of F(q) over the Brillouin 
zone is performed analytically. The same type of method is used to treat the ^2 singularity 
in Eq. fl35l). Here the development of |/3™ n (k, q)| 2 shows that a divergence like ~ also occurs. 
The problem can still be solved by using another function which diverges as - if q — > 0. As 
a matter of fact, this divergence is less severe than the previous one and does not require a 
special treatment since the accuracy of the numerical results is not affected if we neglect this 
divergence. It should be noted that the treatment of the singularity in Eq. (|35|) necessitates 
the evaluation of the symmetrized dielectric matrix for q —>■ 0. As the convergence of the 
head element of this matrix as a function of the number of k-points is slow, the calculation is 
performed separately. All other Brillouin zone integrations are carried out using 10 special 
k-points, which produces accurate numerical results (see the result section). 

The evaluation of the quasiparticle energies requires the determination of the renormal- 
ization factor defined by Eq. (p0|). The derivative of the self-energy is then evaluated by 
using a finite difference scheme with a step equal to 1 eV. The values of the renormalization 
factors are summarized in Table | for the different materials studied here. These constants 
Z are roughly the same for the electron and hole states. Furthermore Z is similar for all the 
materials considered here. The values indicate that account of dynamical renormalization 
is crucial to get quantitatively correct results. At the same time values of Z close to unity 
show that the quasiparticles are well defined and that the GW approximation is reasonable. 

The last point to be discussed is the evaluation of the matrix elements of the exchange- 
correlation potential which appears in Eq. (0). Here some problem arise due to the non- 
linearity of the exchange-correlation potential. Indeed, in the PAW calculation, the exchange 
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correlation potential is calculated by taking into account the valence n v and core electron 
density n c since V x L c DA (r) = V xc DA [n v {r) + n c (r)]. Therefore, the self-energy in the GW ap- 
proximation describe only the valence electrons, and we are obliged to make the assumption 
that the core- valence exchange and core-polarization contributions to the energy of a valence 
state is given by@ 

(^kra | Mrc core— vol |*kn) = (^kn\V xc [n v + n c }\^ kn ) - (#kn|K C H|*kn) (45) 

Such a procedure is some what hidden when pseudo-potential are used since such an oper- 
ation is performed in the unscreening of the pseudo-potential by subtracting 14 c [n„]. The 
shortcoming of this approach is that the ionic pseudo-potential is dependent on the valence 
configuration, reducing the transferability of the potential. Moreover, it has been shown that 
including core corrections to the exchange and correlation is necessary to describe properly 
the structural properties of solidfS. The PAW does not suffer this shortcoming but it seems 
that justifying Eq. fl4"5|) is not an easy task. Considering that the argument mentioned 
before is correct, the quantity which must be subtracted is then defined by 



(^n\V x L c DA [n v (r)}\^ n ) 



^n\V x L c DA [h v (v)}\^ kn ) + EiJ,a<*nk|p?) 

(^v x L c DA K(r)]m - m\v x \ DA [hi{r)]m\ (^ nk > 



(46) 



Table |TTT| shows the matrix elements of the exchange-correlation potential of Si calculated us- 
ing the PAW formalism compared to the results obtained by means of the LAPW method^. 
The agreement between our results and the results of Hamada et al. is remarkably good if 
we consider that these results are based on different calculation schemes and that different 
parameterizations of the exchange-correlation energy are used. 



4- Use of symmetry to reduce computational cost 

If \I/kn is a Bloch wave function solution of the Kohn-Sham equation and R a symmetry 
operation belonging to the point group of the crystal which we suppose to be symmorphic, 
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we can write the way the Bloch wave function transform under such a symmetry operation 

tfknCirt) = J2 D ( R )nAm(r), (47) 
m 

where D(R) nm denote the unitary transformation associated to the symmetry operation R. 
If the state \l/kn is non degenerated, the transformation rule of the wave function simplify 
greatly since D(R) nm = 5 nm . We suppose now that the states considered here are non 
degenerated to simplify the discussion. Using such a relation it can be shown that the 
matrix elements defined by Eq. ( pl| ) satisfy the following relation 

M G m (k, Rq) = M£ m lG (k, q) for Rk = k, (48) 

that means for the symmetry operations belonging to the little group Gk of the point group 
G. Then, it can be proved that the integrand appearing in the Hartree-Fock contribu- 
tion defined by Eq. ([54] ) is invariant under symmetry operations belonging to G^. Such a 
symmetry property reduces the number of q-points for which the integrand has to be calcu- 
lated. Indeed, if BZ^ denotes the irreducible Brillouin zone defined by Gk, the Hartree-Fock 
contribution can be rewritten as 

qeBZ k m occ G Imt^| 

where u>(q) denotes the weight of the q-point. It should be noted here that the term we 
have subtracted from this expression to cancel the Coulomb singularity does not give rise 
to complications since this term is invariant under all symmetry operations of G. Moreover, 
if the state \l/kn is degenerated, we should sum the matrix elements over all degenerated 
states to get a true invariant integrand. The same type of symmetry reduction holds for the 
calculation of the symmetrized static dielectric matrix using the fact that 

M G m ( J Rk,q) = M" m 1 G(k,q) for i?q = q (50) 

and 

e n (Rk) = e n (k) (51) 



We then get 

e GG ,(q,u; - UJ - d GG , - , ,, . 2^ 2^ 2^ 77k — n) — TTk) ^ ' 

The relationship between the matrix elements of e GG / also reduces the computational cost. 
Using the symmetry property of the symmetrized dielectric function e(r, r') = e(Rr, Rr'), it 
can be shown that 

e GG '(Rq,u = 0) = e R -i GR -i G/ (q, to = 0). (53) 

So both the hermiticity of e and the relationship between the matrix elements which results 
from the symmetry operations leaving q invariant are used to reduce the number of matrix 
elements to be computed. Now, we have to remember that the plasmon pole parameters are 
obtained by solving an eigenvalue problem 

J dr' e(r, r') p (q, r') = A p (q) p (q, r) (54) 

By analogy with the resolution of the Schrodinger type equation in a crystal, it can be shown 
that 



(Rq) = <^-i G (q) and A p (i?q) = A p (q) (55) 



These symmetry properties with Eqs. fl29| , pOj, and ^) can be used to show that 



Zp(Rq) = z p (q) and uj p (Rq) = u p (q) (56) 

If the point group G of the crystal does not contain the inversion symmetry, the time reversal 
symmetry could also be implemented. Because of these symmetry relations, the eigenvectors 
and eigenvalues of the symmetrized dielectric matrix are only computed for irreducible q- 
points with respect to the point group of the crystal. Now it can be demonstrated that 
the integrand appearing in the dependent energy contribution to the matrix element of the 
self-energy defined by Eq. ( |3"oD is invariant under symmetry operations belonging to the 
little group of k denoted as in the case of the Hartree-Fock contribution. 
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5. Parallelization 



Because GWA is computationally involved the calculation of QP properties in a se- 
quential computer is time consuming. We have used the message passing interface (MPI) 
to parallelize our numerical code on an IBM SP2. The most straitforward parallelization 
which we were able to perform was any loop or a summation involving k-point mesh over the 
Brillouin zone. First, the irreducible k-point with respect to the point group of the crystal 
were distributed on different processors. Then each processor diagonalize the Hamiltonian 
for a certain number of k-points. It then evaluates the symmetrized dielectric function and 
calculates the ingredients of the plasmon-pole model. Second, each processor calculates the 
self-energy correction for a certain number of k-points. and the results are gathered by 
the root processor to get the GW correction. This simple parallelization scheme made our 
computer code run much fast on the IBM SP2. Our choice of the the k-point mesh was also 
motivated by a future parallelization of the LDA eigenvalue problem for systems with many 
atoms per unit cell. We can set different pools of processors, and each pool will produce the 
diagonalization of a large Hamiltonian in parallel making the code highly scalable with the 
number of atoms of the system. 

III. GW QUASIPARTICLE RESULTS 

A. Quasiparticle results for small and medium band-gap semiconductors: Si, GaAs, 

AlAs, InP, and Mg 2 Si 

In this subsection we present the electronic structure of several small and medium gap 
semiconductors which are used as tests for the implementation of our all-electron PAW-GWA 
method. As mentioned earlier, we have implemented three different types of plasmon-pole 
models available in the literature,01ll@ and compared the resulting quasiparticle energies 
of Si for high symmetry points. Table [Ij shows that these quasiparticle energies are not 
sensitive to the type of plasmon-pole model used to mimic the frequency dependence of the 
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screened interaction. This trend has been confirmed for the other semiconductors studied 
in this paper and was already mentioned in the seminal work of HedinEI on the Jellium 
model. Then we made a detailed comparison with the only full-potential GW calculation 
of Hamada et a/.El where LAPW-GWA calculation of the band structure of Si is presented. 
Table JTJ compares different key ingredients necessary to evaluate the quasiparticle energies 
of Si in the GW approximation. Our calculated matrix elements of the exchange-correlation 
potential (# kn |V^ DA [n^(r)]|\I> kn ) are compared to those obtained by LAPW-GWA. The 
agreement with LAPW is excellent and reflects the accuracy our LDA results obtained by 
PAW method. We have also compared the screened exchange (SEX) and Coulomb-hole 
(COH) contributions to the self-energy with the results of Hamada et a/.El and of Hybertsen 
and Louie!. To make this comparison reliable, we have used the same type of plasmon-pole 
model than Hamada et al. We have found that the agreement of ours results with these 
of Hamada et al. is not fully satisfactorily, but our results are very close to the results 
of Hybertsen-Louie although their calculation is a pseudopotential one and uses another 
type of plasmon-pole model. This observation makes us confident with the PAW-GWA 
results and confirms the fact, as outlined before, that the detailed structure of the screened 
interaction is not important to determine the quasiparticle energies. We should point out 
here that X = YIsex + ^coh which enters the calculation of the quasiparticle energies are 
close to each other whether we consider the PAW, the LAPW, or the Pseudo-Potential 
implementations. 

To calculate the quasiparticle energies it is important to correctly determine the quasi- 
particle renormalization factor Z nk (see Eq. |). As stated earlier, Table | presents Z nk 
calculated for the top valence state at the T point and for the lowest conduction state of all 
semiconductors studied in this paper. These values are in good agreement with the results 
of Hybertsen and Louiei and they seem to be material and state independent and are at the 
vicinity of 0.8. 

Tables [V| § [VII], and |VTTT] represent the calculated PAW-LDA and PAW-GWA 



band energies for the high symmetry points T, X, and L for small and medium band gap 
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semiconductors: Si, GaAs, AlAs, InP, and Mg 2 Si, respectively (The energy scale is relative 
to the top of the valence state maximum). These results are compared to other GWA cal- 
culations obtained using LAPW (for Si) or PP for the other systems as well as to available 
experimental data&^. In these tables we have presented our results for 2 Chadi-Cohen 
k-points0 which correspond to 32 k-points in the whole Brillouin zone (BZ) (results shown 
between parentheses) as well as the most converged values using 10 special k-points^l cor- 
responding to 256 k-points in the BZ. We notice that most of the data in the literature 
are produced with about 32 k-points in the BZ, and are in excellent agreement with our 
unconverged values, however only the converged values should be compared to experiment. 
In our case we have found that the QP eigenvalues are converged with a k-point mesh be- 
yond 256 points. The discrepancy between our GW values and others is often traced to 
differences between LDA values. It's also worth mentioning that Shirley et a/.'s GW results 
for GaAsS agree well with our results (see Table [V]) despite that they are of PP type, and 
that their calculation shows that only the inclusion of core-polarization effect produce gaps 
in agreement with experiment. 

Fig. |3] presents the corresponding band structure along the LT and TX high symmetry 
directions calculated within the LDA and GWA. We notice an overall improvement of the 
excited states eigenvalues compared to these obtained in LDA, whereas the LDA valence 
states eigenvalues are already in good agreement with experiment and GWA results do not 
change this agreement. In all these small and medium band gap semiconductors we remark 
also an average energy shift of the conduction states towards the high energies compared 
to LDA. This energy shift is about the same for Si and GaAs and is about 0.6 eV, and 
increases to about 0.8 eV for AlAs and InP. To be more specific we studied the range of 
applicability of the so called scissors-operator shift. To this end we have evaluated the 
deviation of the difference of LDA and GWA direct band gap determined at the T, X, 
and L symmetry points. We have found that these deviation of the GWA and LDA energy 
differences are the lowest for Si, and Mg 2 Si. The maximum deviation is about 0.06, and 0.04 
eV, respectively, and it occurs from L to V in both materials. The deviations are somewhat 
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larger for GaAs, AlAs, and InP, and the maximum deviation is about, 0.16, 0.15, and 0.13 
eV, respectively. They are all from T to X. These small deviations indicate that the GWA 
does not change much the LDA dispersion across the Brillouin zone, justifying the use of 
the scissors-operator shift for the calculation of the dielectric function for small and medium 
band gap semiconductorsSHfl. We will see in the next subsection that such deviations are 
much larger in absolute values for wide-band-gap semiconductors. 

As for Mg2Si we believe that it is the first time that this compound is studied within the 
GW approximation. The PAW-LDA and PAW-GWA band energies for the high symmetry 
points T, X, and L are shown in table [VIII| . Due to the lack of photoemission experiments, 
the GW results are compared with optical measurements, making the assumption that exci- 
tonic effects are negligible. The GWA results are in good agreement with the experimental 
results and compare favorably with the EPM calculation of Au-Yang et al.0 Fig. ^ presents 
the corresponding band structure along the LT and TX high symmetry directions within 
the LDA and GWA. 

Fig. |] shows the LDA and GWA calculated minimum band gaps for all studied semicon- 
ductors, and are compared to the experimental results. A perfect agreement with experiment 
is achieved when the calculated value falls on the dashed line. We notice that for most of 
the small and medium gap semiconductors GWA does not account for the whole correction 
of the band gap. The disagreement with experimental band gaps is most probably due to 
the procedure used to decouple the core and valence electrons. Indeed, we used the LDA 
exchange-correlation potential with the valence electron density to estimate the LDA coun- 
terpart of the self-energy and we believe that a more refined treatment should rely on the 
evaluation of the core- valence exchange interaction within the Hartree-Fock approximationS 
However, this latter suggested procedure is difficult to implement in our formalism and the 
outcome would not necessarily improve upon the former simple scheme. On the other hand, 
it is interesting to mention that a first order vertex and self-consistent corrections to the 
RPA polarizability and to the self-energy within the GWA increase the direct energy band 
gaps of Si at the T, L, and X points by about 0.36, 0.44, and 0.39 eV, respectivelyS. It 
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seems then that there is no compensation between the vertex correction and the selfconsis- 
tency as original thought. If we start from our converged Si GWA results, and use these 
latter results we could improve the agreement between our calculated band gaps and ex- 
perimentent. However, when some PP-GWA results are used, the additional corrections 
overestimates the experimental gaps. At the present time it looks like the question of the 
band gaps is not fully solved. It is also important to remark that these vertex corrections 
and the corrections arising from the selfconsistency are not quite accurate since the starting 
point of these calculations is a noninteracting Green's function instead of the selfconsistent 
Green's function as suggested by Hedini. 



B. Quasiparticle results for wide band gap semiconductors and insulators: C, and 

LiCl 

It is of interest to compare all-electron GWA calculations for wide band gap semiconduc- 
tors and insulators to existing PP calculations. Wide band gap semiconductors are somehow 
puzzling in contrast to small and medium semiconductors: While the LDA band gap of these 
materials are significantly underestimated compared to experiment, the LDA static dielectric 
function are usually in good agreement with experimental results, see for instant Ref.0@. 

Tables [EX|, and |X| show the calculated PAW-LDA and PAW-GWA band energies for 
the high symmetry points T, X, and L for wide band gap semiconductors: C, and LiCl, 
respectively (The energy scale is relative to the top of the valence state maximum). These 
results are compared to other GWA calculations obtained using PP-GWA methodic and 
to experimental data whenever availableH'&i^l. Fig. [| presents the corresponding band 
structure along the LT and TX high symmetry directions within the LDA and GWA. The 
C and LiCl values have been computed using 32 k-points as well as 256 k-points in the 
Brillouin zone. The results of the two sets of k-points are in good agreement showing that 
the set of 32 k-points is good enough. For C the calculated QP eigenvalues are in good 
agreement with experiment and the PP calculations. For LiCl only the experimental band 
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gap is available and is slightly larger than our GW value. It's worth mentioning that we 
didn't update the Green function to get our GW values. Such a procedure leads to an 
increase of the GW band gap by about 0.3 eV and then to a better agreement between our 
results and PP results of Hybertsen and Louie!. For these wide-gap materials we looked 
also to the applicability of the scissors-operator shift. We calculated the maximum change 
of the difference between the GWA and LDA direct band gap at the T, X, and L symmetry 
points. We found that the maximum deviation across the Brillouin zone is for C and is 
about 0.32 eV. It occurs for the L and T differences, whereas it is about 0.28 eV for LiCl 
and occurs for X and T energy differences. These deviations seem to be somewhat larger 
in absolute values (about twice the value found for GaAs) than for small and medium-gap 
semiconductors, and may make the use of the scissors-operator shift, for the computation 
of the optical properties, less applicable. However, if we compare these energy deviations to 
the size of the band gap, we find that the largest ratio occurs for C and is only 5%, whereas 
it is about 10% for a medium gap semiconductor such as GaAs. 



IV. CONCLUSION 

We have implemented a GWA based on an all-electron method using the recently devel- 
oped projector augmented wave (PAW) method.@ The knowledge of the one-electron Green 
function of the PAW Hamiltonian allows us to construct the quasiparticle self-energy within 
the GWA, in which the dynamical screening of the electron-electron interaction arises from 
a plasmon model dielectric function0Elil for which the parameters are adjusted to the di- 
electric function calculated using the random-phase-approximation (RPA). We have tried 
various plasmon model dielectric functions for the screening of the Coulomb interaction and 
showed that the quasiparticle energies are insensitive to the type of the model used. 

Using this new GWA method, we have determined for the first time the GWA quasi- 
particle electronic structure of Mg2Si. We have found that our LDA results are in good 
agreement with the empirical pseudopotential, and that the GWA shifts almost rigidly the 
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empty states by about 0.32 eV towards higher energies. 

Concerning the other semiconductors studied here, we have found an overall agree- 
ment of our calculated electronic structure of various semiconductors with existing GWA 
pseudopotential calculations performed by different groups@iE00, and with experimental 
result However, for detailed comparisons, the all-electron GWA band 
gaps are slightly smaller that most of the PP results. Nevertheless, some of the PP results^ 
are much closer to ours. One of the possible source of the the all-electron results compared 
to the PP ones is the way the decoupling of the core and valence electrons is done. In the 
PP approach this decoupling is some what hidden. In fact, we have used the LDA exchange- 
correlation potential with the valence electron density to remove the LDA counterpart of 
the self-energy. We believe that a more refined treatment should rely on the evaluation of 
the core- valence exchange interaction within the Hartree-Fock approximationi which we find 
difficult to implement in our current all-electon PAW-GWA approach. On the other hand, 
the first order vertex and self-consistent corrections to the RPA polarizability and to the 
self-energy within the GWA are shown to increase the direct energy band gaps of Si at the 
T, L, and X points by few tenths of an It seems then that there is no compensation 

between the vertex correction and the selfconsistency as it was always assumed. However, 
these corrections are not quite accurate since the starting point of these calculations is a 
noninteracting Green's function instead of the selfconsistent Green's function as suggested 
by Hedini. 

To our knowledge this is the first all-electron GWA calculation that has corrected LDA 
eigenvalues for three type of semiconductors: small, medium, and wide band gap, and that 
questioned the correctness of the band gap of semiconductors obtained by means of PP- 
GWA. We hope that this work would be used as a reference and triggers off further interest 
on an all-electron GW approach. 
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TABLES 

TABLE I. The renormalization constants Z for the hole state at the top of the valence band 

(VBM) and the electron state near the bottom of the conduction band (CBM) for C, Si, GaAs, 
AlAs, InP and SiMg 2 . 

C Si GaAs AlAs InP Mg 2 Si 

Z \bm °- 85 °- 80 °- 80 °- 81 °- 79 °- 76 

Z Cbm °- 87 °- 81 °- 81 °- 82 °- 82 °- 79 
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TABLE II. Quasiparticle (QP) energies of Si for several states (in eV) and for three different 
types of plasmon-pole models. We notice that the QP energies are less sensitive to the type of 
plasmon pole model used. 



Plasmon-pole model 
Von der Linden and Horsch a Hamada et aP 



Engel-Farid c 



P25'u 
Tl5c 



-12.01 
0.00 
3.09 
3.88 



-11.99 
0.00 
3.07 
3.86 



-11.89 
0.00 
3.12 
3.86 



-7.97 
-3.00 
1.01 



-7.96 
-3.00 
0.99 



-8.01 
-3.01 
1.04 



Lie 
Lie 



-9.73 
-7.21 
-1.25 
1.96 
3.85 



-9.72 
-7.20 
-1.25 
1.94 
3.83 



-9.75 
-7.19 
-1.29 
1.94 
3.86 



a Ref.§, b RefM c R e f.@ 
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TABLE III. Matrix elements of the exchange-correlation potential and of the screened exchange 
(SEX) and the Coulomb hole (COH) contributions to the self-energy of Si for several states (in 
eV). The Vxc is i n excellent agreement with the LAPW result Jll, whereas the SEX and COH are 
much closer to the PP work of Hybertsen and Louie! and disagree with the LAPW results. 
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TABLE IV. Quasiparticle energies of Si for several states (in eV). The calculation of the 
self-energy is performed using 10 special k-points, 200 bands and 283 reciprocal lattice vectors. 
The size of the polar izability matrix is 137x137 and the plasmon pole model of von der Linden 



and Horscl 


£3 is used. 
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TABLE V. Quasiparticle energies of GaAs for several states (in eV). The calculation of the 
self-energy is performed using 10 special k-points, 200 bands and 307 reciprocal lattice vectors. 
The size of the polar izability matrix is 169x169 and the plasmon pole model of von der Linden 
and Horsch0 is used. Our results are compared with those of Rohlfing et al0 and those of Shirley 
et al§. The results of Shirley are given without (with) core polarization effects. Here E g is the 
minimum band gap. 
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TABLE VI. Quasiparticle energies of AlAs for several states (in eV). The calculation of the 
self-energy is performed using 10 special k-points, 200 bands and 283 reciprocal lattice vectors. 
The size of the polar izability matrix is 169x169 and the plasmon pole model of von der Linden 
and Horsch0 is used. Our results are compared with those of Godby et a!0 (GS) obtained using 
the RPA dielectric function for the screening of the Coulomb interaction, and these of Shirley, Zhu, 
and Louie (SZL) using a plsmon-pole model. The Godby et al results include spin-orbit coupling 
and the lower energy of a spin orbit pair is shown between brackets. Here E g is the minimum band 



gap. 
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TABLE VII. Quasiparticle energies of InP for several states (in eV). The calculation of the 



self-energy is 


performed usin; 


I 10 special k-points, 200 bands and 331 


reciprocal lattice vectors. 


The size of the polarizability matrix is 


169x169 and the plasmon pole 


model of 


von der Linden 


ITT 1 h (1 

and Horscntil 


is used. Our results are 


compared with those of Hott0. 
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TABLE VIII. Quasiparticle energies of Mg2Si for several states (in eV). The calculation of the 
self-energy is performed using 2 special k-points, 200 bands and 645 reciprocal lattice vectors. The 
size of the polarizability matrix is 113x113 and the plasmon pole model of von der Linden and 
Horsch0 is used. Here E„ is the minimum band gap. 
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TABLE IX. Quasiparticle energies of C for several states (in eV). The calculation of the 
self-energy is performed using 10 special k-points, 200 bands and 387 reciprocal lattice vectors. 
The size of the polar izability matrix is 169x169 and the plasmon pole model of von der Linden 
and Horsch0 is used. Here E g is the minimum band gap. 
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TABLE X. Quasiparticle energies of LiCl for several states (in eV). The calculation of the 
self-energy is performed using 10 special k-points, 200 bands and 331 reciprocal lattice vectors. 
The size of the polar izability matrix is 259x259 and the plasmon pole model of von der Linden 
and Horsch0 is used. Here E g is the minimum band gap. 
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FIGURES 




wave vector 



FIG. 1. Calculated Engel-Farid Plasmon model band structure of Si along L, T, and X 
high symmetry directions. The agreement with the results of Engel and FaridS, and AulbuiS 
are excellent. For small k wave vectors the lowest plasmon band shows a quadratic dispersion 
(up-triangle curve) u>o(k) = ujq(0) + o;|k| 2 , with a dimensionless direction-dependent dispersion 
coefficient a. We find cjo(O) = 15.7 eV and «a = 0.33 in good agreement with the values of 15.91 
eV and 0.34 of Engel and Farid0 as well as the experimental values^! of 16.7 eV and 0.41. 
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FIG. 2. Ab-initio calculated real and imaginary parts of the inverse dielectric function (dashed 
curve) of Si compared to the plasmon pole model of Hamada et al (solid line) and with experiment 
(open circles)@. 
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FIG. 3. Calculated electronic band structures along high symmetry directions for some small 
and medium band gap semiconductors: Si, GaAs, AlAs, InP, and SiMg2 (in eV). The dashed lines 
display LDA results calculated with an energy cut-off of 15Ry (cf. Table fV|,|V|JV^JVI^ ; and |vn| 
column 2). The solid lines with dots show the GW results based on these LDA results (cf. Tables 
[Tv|,|V|, [V|,[Vn|, and [Vll| column 3). The energy scale is relative to the top of the valence state 
maximum (VBM). 
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FIG. 4. Calculated LDA and GWA band gap compared to experimental results. The filled 
circles represent the LDA values, the open circles the GW values with 10 special k-points (cor- 
responding to 256 points in the BZ) and the up-triangles the GW values using 2 chadi-Cohen 
k-points (corresponding to 32 points in the BZ). A perfect agreement with experiment is achieved 
when the calculated value is on the dashed line. 
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FIG. 5. Calculated electronic band structures along lines of high symmetry for some large 
band gap semiconductors: C, and LiCl (in eV). The dashed lines display LDA results calculated 
with an energy cut-off of 45Ry (cf. Table |x||x|, column 2). The solid lines with dots show the GW 
results based on these LDA results (cf . Tables |X|,[X|, column 3) . The energy scale is relative to the 
top of the valence state maximum (VBM). 
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